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Abstract. In a recent paper, Edvardsson et a/. (2002) propose a new solution for the spin evolution of the Earth 
and Mars. Their results differ significantly with respect to previous studies, as they found a large contribution 
on the precession of the planet axis from the tidal effects of Phobos and Deimos. In fact, this probably results 
from the omission by the authors of the torques exerted on the satellites orbits by the planet's equatorial bulge, 
as otherwise the average torque exerted by the satellites on the planet is null. 
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1. Introduction 

In a recent paper, Edvardsson et aL(2002) propose a new 
solution for the spin evolution of the Earth and Mars. 
Considering the absence of a precise evaluation of the er- 
rors due to the integrator and the absence of relativity in 
the model, one could discuss the use of the word 'accurate' 
in the title of the paper, but unfortunately there are some 
more important flaws in this paper. In their integration of 
the spin of Mars, the authors found that the integration 
of the spin changes in a large amount when Phobos and 
Deimos are taken into account (see their Fig. 12). They also 
notice that their " curve without the moons is very similar 
to the curve given by Bouquillon & Souchay (1999) " , who 
included the moons. In fact, the paper of Edvardsson et 
a/.(2002) is largely in error on this point. 



2. Precession due to a distant satellite 

Let us consider a planet V with momentum of inertia 
A = B < C orbiting the Sun on a fixed ellipse (we will not 
consider here the planetary perturbations or the pertur- 
bations due to the satellite presence), and a satellite S of 
mass m orbiting the planet. Let (i,j,k) be a basis linked 
to V with k associated to the axis of maximum inertia 
C. Let (I,J,K) be a fixed reference frame, with origine 
in the direction of I, and K normal to the orbital plane 
of V (Fig.l). a and e are the semi- major axis and eccen- 
tricity of 5, n and i the longitude of the axcending node 
of the satellite orbits over the orbital plane of the planet 
(I, J), while V is the true longitude, and u the argument 
of perihelion. If r is the radius vector from the planet's to 



Equator{t) 




Fig. 1. 



the satellite's center of mass, with modulus r and unitary 
vector V = r/r, the torque exerteed by S on V is 



-V AXv 



(1) 



where G is the gravitational constant and X the ma- 
trix of inertia (I = diag{A, A,C)). Noting that T = 
dia.g{A,A,A) — dia,g{0,0,C—A), (1) can also be expressed 
as (see Murray, 1983) 



3Gm 

F = 3-(C - A)k A (vv^)k 



(2) 
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where v"^ denotes the transposed of v (v'^k is thus the 
dot product of the two vectors v and k). We will average 
F over the fastest angle of this problem, that is over the 
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mean anomaly M of the satellite (the rotational angle 9 
is removed after the asumption A = B). As k does not 

depends on A^f, the only expression to average in (2) is 
7 = w'^/r^. In the (I, J,K) frame, the coordinates of v 
are 



(cos(w + w) ' 
sin('u + u) 

with 

/I 

TZi (a) = cos a — sin a | ; 
\ sin a cos a 

cos a — sin a ^ 
7?,3(a) = I sin a cos a 
1 



(3) 



(4) 



As dv = a^Vl — e^dM, after averaging over M, we 
obtain 



< 7 >M= 



1 



:[Id-M] 



2a3(l- 6^)3/2 
where M. is expressed in (I, J, K) as the 3x3 matrix 



(5) 





= 7^3(^^)7^l(^) I o o 1 7^l(-^)7^3(-^^) (6) 

1^ 

When the precession of the node is rapid with respect 

to the procession of the spin axis of the planet (as for the 
Moon around the Earth), one can also average over fi. We 
have thus in the basis (I, J, K) 

^ \0 1/ 

Finally, as the parts involving the identity Id cancel, 
we are left with 



^GmiC-A) , 3sin2i 
< ' = - 2a3(l -e^)3/- (^ - -^^"^ ' 

- ^^-(^-^W-^)cosskAK 



(8) 



2a\l - e2)3/2 



which leads to the classical contribution of the preces- 
sion of spin axis k around the normal to the orbital plane 
K, when the rotational rate v is large with respect to the 
precession frequency p. 

3Gm C-A^ Ssin^i 
= -2.a3(l-eT/^ C 

3. Precession due to a close satellite 

The previous study is in fact valid uniquely in the case 
of a satellite sufficiently far from the planet, so that the 
precession of its orbit is mostly driven by the solar pertur- 
bations, and not by the torque exerted by the equatorial 



bulge of the planet. In the case of a close satellite, it was 
shown (Goldreich, 1965, Kinoshita, 1993) that instead of 
precessing around the normal K to the orbital plane (I, J), 
with a roughly constant inclination i with respect to (I, J), 
the satellite will precess with respect to k with a nearly 
constant inclination i' with respect to the equator plane 
(i, j). In this case, all the previous study is still valid, but 
the orbit of the satellite is now referred to the equatorial 
reference frame (i, j, k), and the second averaging is made 
with respect to the longitude of the node, ft', of the satel- 
lite orbit with respect to the equatorial plane (i, j). The 
average expression (6, 7) is still valid, when replacing Q 
and i by f2' and i'. The average < M >a' will then be 
given by (8), after changing i with i', but this expression is 
now obtained in the reference frame of the planet (i, j, k). 
As k A kk'^k = k A k = 0, we have now for the averaged 
torque 



< r >M,n'= 0. 



4. Discussion 



(10) 



Thus, contrarily to what is found by Edvardsson et 
o/.(2002), the averaged torque exerted by Phobos and 
Deimos on Mars should be null (at first order). The ori- 
gin of the authors' error can be traced by looking to their 
figures 1 and 2. The inclination of the satellites with re- 
spect to the equator start with a value close to zero at 
the origin, and then increase up to about 50 degrees and 
then back to zero with a period of about 6000 and 1500 
years for Phobos and Deimos, while a back on the en- 
veloppe calculation gives respectively 5400 and 1400 years 
for the precession of the node of these satellites submitted 
uniquely to solar perturbation. The issue is then clear. By 
neglecting the torque exerted on the satellites orbits by 
the equatorial bidge of Mars, the authors found that in- 
stead of remaining close to the equatorial plane of Mars, 
as demonstrated by Goldreich (1965), the satellites were 
precessing along the orbital plane. Instead of finding a 
zero value for the mean precession torr^ue (11), they ob- 
tain an efi'ect following Eq. (10) which gives a contribution 
of about 0.5" /year to the precession of mars axis, in agree- 
ment with the difference observed by the authors in figure 
12. It is surprising that the authors did not wonder how 
Phobos and Deimos could both be at present with an in- 
clination of less than 2 degrees with respect to the equator 
if they were precessing along the orbit. Goldreich indeed 
realised more that this 'would amount to an unbelievable 
coincidence'. 
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